Accepted for Publication in ApJ 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



00 

O 
O 

(N 

C 

a 
>—> 

OS 



6 



> 
(N 



O 

oo 
O 



A PRECISE PROPER MOTION FOR THE CRAB PULSAR, AND THE DIFFICULTY OF TESTING 
SPIN-KICK ALIGNMENT FOR YOUNG NEUTRON STARS 

D. L. Kaplan 1 , S. Chatterjee 2 , B. M. Gaensler 2 ' 3 , and J. Anderson 4 

Accepted for Publication in ApJ 

ABSTRACT 

We present a detailed analysis of archival Hubble Space Telescope data that we use to measure the 
proper motion of the Crab pulsar, with the primary goal of comparing the direction of its proper 
motion with the projected axis of its pulsar wind nebula (the projected spin axis of the pulsar). 
Combining data from 47 observations spanning > 10 yr with two different instruments, and using the 
best available measurement techniques and latest distortion models, we are able to demonstrate that 
our measurement is robust and has an uncertainty of only ±0.4 mas yr" 1 on each component of the 
proper motion. However, we then consider the various uncertainties that arise from the need to correct 
the proper motion that we measure to the local standard of rest at the position of the pulsar and find 
/i Q = — 11.8±0.4±0.5masyr _1 and fis = +4.4±0.4±0.5masyr _1 relative to the pulsar's standard of 
rest, where the two uncertainties are from the measurement and the reference frame, respectively. If 
we then wish to compare this proper motion to the symmetry axis of the pulsar wind nebula, we must 
consider the unknown velocity of the pulsar's progenitor (assumed to be ~ lOkms" 1 ), and hence add 
an additional uncertainty of ±2 masyr -1 to each component of the proper motion, although this could 
be a factor of 10 larger if the pulsar's progenitor had an anomalously high velocity (> 100 krns" 1 ). This 
implies a projected misalignment with the nebular axis of 14° ± 2° ± 9°, consistent with a broad range 
of values including perfect alignment. We use our proper motion to derive an independent estimate for 
the site of the supernova explosion with an accuracy that is 2-3 times better than previous estimates. 
We conclude that the precision of individual measurements which compare the direction of motion of 
a neutron star to a fixed axis will often be limited by fundamental uncertainties regarding reference 
frames and progenitor properties. The question of spin-kick (mis)alignment, and its implications for 
asymmetries and other processes during supernova core-collapse, is best approached by considering a 
statistical ensemble of such measurements, rather than detailed studies of individual sources. 
Subject headings: astrometry — pulsars: individual (PSR B0531+21, Crab) — stars: neutron 



1. INTRODUCTION 

The high velocity nature of the neutron star 
population has been apparent almost as long as 
the ex istence of neutr on stars has been recog- 
nized ()Gunn fc Ostrikeil Il970fh Recent statisti- 
cal studi es of the radio pulsar velocity distribu- 
tion (e.g.. lArzoumanian et al . 2002: Brisk en et all [2003: 
iHobbs et al.ll2005t iFaucher-Giguere fc Kaspill2006[ ) yield 
mean three-dimensional velocities of 300-500 kms -1 , 
with a high velocity tail extending beyond 1000 kins" 1 . 

A variety of physical mechanisms have been pro- 
posed as the origin of high velocities. Perhaps the 
first was d isruption of b i naries through mass loss in su- 
pernovae (iBlaauwl 119611 : iGott. Gunn. fc Ostrikerl 119701 : 
llben fc Tutukovl 11996ft although it is difficult for bi- 
nary disruption alone to account for some of the 
highe st observed velocitie s (lHarrison. Lyne. fc Anderson! 
119931: IChatteriee et al.l l2005f ). The most natu- 
ral source of such high velocities appears to be 
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asymmetries in the birth supernovae of pulsars 
( Shkl ovskiil [19701: Ivan den Heuvel fc van ParadiiH 119971 : 
iPortegies Zwart fc van den Heuvell Il999f ). although ex- 
actly how an asymmetry in the core collapse process 
in a supernova explosion is converted to a birth kick 
imparted to a nascent neutro n star remains unclear 
(|Lai. Chernoff. fc Corded 120011 ) . While hydrodynamic 
or co nvective instabilities are t he most plausible route 
e.g.. iBurrows fc Havel 119961: IJanka fe Mueller] 119961: 



Lai fc Goldreichl 120001: IScheck et all l2004t IJanka et al l 
20051: IScheck et all 120061 : iBurrows et al.ll2006h . more ex- 



otic mechanisms such as asymmetric neutr ino emission 
in th e presence of strong magnetic field s (lArras fc Lai 
1999h or some combination of the above (jSocrates et al 



2005) cannot be ruled out. 

Of these kick mechanisms, many predict kicks 
vectors that relate to the spin axis orientation of 
the nascent neutron star: the alignment (or lack 
thereof) of the natal kick with the neutron star spin 
axis could provide a specific discriminant between 
various mechanism s (Burrows. Haves, fc FrvxeU H9951 
[ Spruit fc Phirmevl [199^ iCowsikl [19981 : lLai et all 
lRomanill2005f) . Even such parameters such as the num- 
ber and timescale of kick components, coupled with 
the initial spin period of the neutron star, can be 
constrain ed through observations of an ensemble o f 
sources (Deshpande. Ramachandran, & Radhakrishna: 
19991: IJohnston et all 120061: IWang et alJ I2006L 1200 



Ng & Romani 2007; Rankin 2007 



2 



1.1. The Crab Pulsar and its Nebula 

The Crab pulsar (PSR B0531+21) and nebula have 
been observed by virtually every telescope capable of 
pointing at the system, and the Crab may be the most 
studied system in all of astronomy. The Crab nebula 
possesses a general symmetry axis, visible in images at 
most wavelengths. Recent observati ons have delineated 
this axis with strik ing clarity (e.g., [Hester et all 120021 : 
iNg fe Romanl[200l . The X -ray jet, in particular, allows 
us to trace the symmetry axis to the neutron star loca- 
tion, and provides a natural association with the rotation 
axis of the pulsar itself (since every other vector would 
be rotation averaged). The symm etry axis is roughly 
align ed with the proper motion (e.g. JCaraveo fc Mignaiill 
Il999f >. but the observational uncertainties have made the 
alignment hard to quantify. A precise proper motion vec- 
tor for the Crab pulsar could be quantitatively compared 
to the jet direction to establish whether the natal kick is 
aligned with the spin axis, as many theories predict. 

Given its prominent place in our understanding of neu- 
tron stars, it is perhaps surprising that the proper motion 
and distance of the Crab pulsar are not better known. 
Compared to many fainter objects, the precision of our 
measurements is lacking. While there were a number 
of early attempts to measure the proper motion of the 
Crab pulsar, these were generally inconsistent with each 
other ([Minkow ski 1970) . Perhaps the first reli able mea- 
surement was that o f iWvckoff fc Murravl (fl977l . hereafter 
IWM77D . who found 5 (n a ,vs) = (-13±2, 7±3) masyr" 1 
from photographic plates spanning epochs from 1899 to 
1976. There have not been any direct (i.e. geometric) 
distance measurements of the pulsar itself, but the dis- 
tance was estimated based on various lines o f evidence 
to lie between 1.4 and 2.7kpc (| Trim ble 1973). In spite 
of the wealth of observations since then, including a 
treasure trove of Hubble Space Telescope (HST) images 
with a time baseline of > 10 years, these estimates had 
not significantly im p roved until rece ntly. F or example, 
ICaraveo fc Mignanil ([1999I . hereafter ICM99I ) estimate a 
proper motion (/x a , /is) = {— 17±3, 7±3) masyr -1 , which 
is consistent with the earlier estimate, but does not im- 
prove upon its accuracy. 

The main obstacle to more accurate measurements is 
not (as in many cases) the limitations of faint objects, 
but rather the fact that the Crab nebula is too bright 
for most interferometric radio observations (it raises the 
system temperature too much), that there are no suit- 
able interferometric calibrators nearby, and that the ro- 
tational stability is not sufficient (due to glitches) for 
precise pulse timing over a long time baseline. Because 
of these reasons, high angular resolution optical observa- 
tions are so far the only way to measure the astrometric 
parameters (proper motion and parallax) of the Crab 
pulsar, and with the current generation of instruments 
w e are limited to data from HST 

INg fc RomanH (|2006l hereafter iNROll attempted a sig- 
nificantly more det ailed as trometric analysis of the Crab 
pulsar compared to |CM99l taking advantage of new HST 
data spanning 7 years and trying to acc ount fo r many 
sources of uncertainty not addressed by CM99. NR06 

J The analysis of WM77 was based on the B1950 frame, but at 
this level of precision the precession between that and the J2000 
frame does not change the proper motion significantly. 



found a result that is discrepant with that of ICM99I : 
(Mc^O = (-15.0±0.8,1.3±0.8)masyr" 1 . This shows a 
significant misalignment with the projected spin-axis of 
the pulsar (26° ±3°) and as such reverses previously held 
notions of spin -kick alignment, but as we discuss below 
(and as lNR06l acknowledge) even this analysis still is not 
as accurate as possible. Additionally, a large number of 
new observations with the Advanced Camera for Surveys 
(ACS) have become publicly available. These were taken 
primarily for studying the dynamics and polarization o f 
the Crab's pulsar wind nebula fe.g.. iHester et alj '2002). 
and as such they were not ideal for astrometry (the ex- 
posure times were long enough that the pulsar saturated, 
the dithering strategy was not optimal, and they used a 
limited range of roll angle), but nonetheless they are an 
important resource. 

Motivated by the importance of the Crab pulsar in our 
understanding of neutron stars and supernova remnants 
in general, by the large amount of data available on it, 
and by the limitations of previous analyses, we have at- 
tempted to re-measure the proper motion of the Crab 
pulsar as well as assess the possibility of a parallax mea- 
surement with future data (although given the recent 
failure of the ACS instrument, future observations may 
not be possible until the installation of the upcoming 
Wide Field Camera 3). 

The organization of our paper is as follows: in 
§ [5] we describe our analysis, noting departures from 
previous analyses, although the majority of the fit- 
ting techniques are similar to those we used in 
iKaplan. van Kerkwiik. fc Anderson! (|2007f ). and we refer 
readers there for more details. After refining the proper 
motion measurement, we discuss in ? 12.51 the limitations 
on our knowledge of the proper motion imposed by the 
unknown velocity of the pulsar's progenitor, as well as 
uncertainties in the corrections to the pulsar's local stan- 
dard of rest. These transformations and their associated 
uncertainties limit the accuracy of the comparison be- 
tween the proper motion and the projected spin-axis of 
the pulsar. We then give our conclusions in § [3J Finally, 
we include a discussion of the prospects for a parallax 
distance for the Crab pulsar in Appendix [AJ In what 
follows, we define our proper motions in Right Ascen- 
sion and Declination (fi a ,fis) so that the scales are the 
same and no cos S term is necessary. All uncertainties 
are 1-a unless otherwise stated, and all position angles 
are measured east of north. 

2. ANALYSIS 

We started our analysis by examining the available 
archival HST data for the Crab pulsar. Twenty ob- 
servations using the Wide Field and Planetary Cam- 
era 2 (WFPC2) with the F547M filter (a filter centered 
at y~band but somewhat narrower, designed to avoid 
bright emission lines) spanning ~2 years were analyzed 
by[CM99, who were able to determine a proper motion 
for the pulsar th at agreed with that obtained from the 
ground (|WM77f ). However, the precisions of both of 
those measurements were limited and we have a num- 
ber of reasons to suspect the analysis of the HST data. 

When we examined the HST data used in the prior 
analyses in detail we noticed that the pulsar itself was 
very saturated. This is not unexpected: using the 
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TABLE 1 

Observation Summary 



Pair Root a MJD Date Instrument/ Exp. Crab (x,y) ra „ b PA N sta .r B c NR06 

Number Name Detector Group d 

(sec) (pixels) (deg.) 



1 


u2bx05 


49723.7 


1995-Jan-07 


WFPC2/WF3 


800.0 


162.98 


130.19 


309.0 


11 


1 


2 


u2bx05 


49723.8 


1995-Jan-07 


WFPC2/WF3 


1000.0 


150.84 


117.81 


309.0 


11 


1 


3 


u2u601 


49943.7 


1995-Aug-15 


WFPC2/PC 


1000.0 


470.47 


371.54 


-48.6 


5 




4 


u2u602 


50026.6 


1995-Nov-06 


WFPC2/PC 


1000.0 


340.13 


300.08 


-25.3 


1 




5 


u2u603 


50080.4 


1995-Dec-29 


WFPC2/PC 


1000.0 


296.63 


600.35 


128.7 


1 


2 


6 


u2u604 


50102.3 


1996-Jan-20 


WFPC2/PC 


1000.0 


297.29 


600.99 


128.7 


4 


2 


7 


u2u605 


50108.3 


1996-Jan-26 


WFPC2/PC 


1000.0 


297.06 


603.05 


128.7 


4 


2 


8 


u2u606 


50114.5 


1996-Feb-02 


WFPC2/PC 


1000.0 


297.83 


600.07 


128.7 


1 


2 


9 


u2u607 


50135.3 


1996-Feb-22 


WFPC2/PC 


1000.0 


297.64 


601.39 


128.7 


4 


2 


10 


u2u608 


50189.6 


1996-Apr-17 


WFPC2/PC 


1000.0 


265.89 


598.65 


128.7 


4 


2 


11 


u61m01 


51580.1 


2000-Feb-06 


WFPC2/WF3 


1100.0 


372.75 


273.44 


312.7 


17 


3 


12 


u61m02 


51589.9 


2000-Feb-16 


WFPC2/WF3 


1100.0 


372.77 


273.30 


312.7 


19 


3 


13 


u61m03 


51600.2 


2000-Feb-26 


WFPC2/WF3 


1100.0 


371.49 


270.99 


312.7 


19 


3 


1 1 


u61m04 


51610.5 


2000-Mar-08 


WFPC2/WF3 


1100.0 


372.25 


272.79 


312.7 


17 


3 


15 


u61m05 


51620.4 


2000-Mar-17 


WFPC2/WF3 


1100.0 


372.66 


273.40 


312.7 


19 


3 


16 


u50v04 


51796.9 


2000-Sep-10 


WFPC2/WF3 


918.0 


216.87 


342.64 


132.7 


8 


5 


17 


u50v05 


51809.0 


2000-Sep-22 


WFPC2/WF3 


1090.0 


216.78 


342.43 


132.7 


9 


5 


18 


u50v06 


51818.8 


2000-Oct-02 


WFPC2/WF3 


872.0 


216.88 


342.45 


132.7 


8 


5 


19 


u50v07 


51829.9 


2000-Oct-13 


WFPC2/WF3 


1200.0 


217.17 


341.14 


132.7 


9 


5 


20 


u50v08 


51840.8 


2000-Oct-24 


WFPC2/WF3 


916.0 


216.75 


342.55 


132.7 


7 


5 


21 


u50vl0 


51863.0 


2000-Nov-15 


WFPC2/WF3 


889.5 


216.78 


342.48 


132.7 


7 


5 


22 


u50vll 


51873.1 


2000-Nov-25 


WFPC2/WF3 


1100.0 


225.62 


334.16 


132.7 


8 


5 


23 


u50vl2 


51884.4 


2000-Dec-06 


WFPC2/WF3 


1000.0 


216.78 


340.60 


132.7 


8 


5 


21 


u50vl3 


51896.3 


2000-Dec-18 


WFPC2/WF3 


1000.0 


371.97 


274.06 


312.7 


17 


6 


25 


u50vl4 


51906.7 


2000-Dec-29 


WFPC2/WF3 


1000.0 


371.88 


273.82 


312.7 


19 


6 


26 


u50vl5 


51918.5 


2001-Jan-09 


WFPC2/WF3 


1200.0 


387.46 


275.98 


312.7 


19 


6 


27 


u50vl7 


51939.4 


2001-Jan-30 


WFPC2/WF3 


1200.0 


387.68 


276.10 


312.7 


18 


6 


28 


u50vl8 


51950.5 


2001-Feb-10 


WFPC2/WF3 


1200.0 


383.72 


278.22 


312.7 


17 


6 


29 


u50vl9 


51961.5 


2001-Feb-21 


WFPC2/WF3 


1200.0 


388.42 


276.09 


312.7 


17 


6 


30 


u50v20 


51972.8 


2001-Mar-05 


WFPC2/WF3 


1200.0 


389.02 


275.64 


312.7 


18 


6 
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32 


u50v22 


51994.5 


2001-Mar-27 


WFPC2/WF3 


1200.0 


389.81 


275.43 


312.7 


18 


6 


33 


u50v23 


52005.2 


2001-Apr-06 


WFPC2/WF3 


1200.0 


389.98 


275.55 


312.7 


17 


6 


34 


u50v24 


52016.3 


2001-Apr-17 


WFPC2/WF3 


1200.0 


389.50 


275.67 


312.7 


18 


6 


35 


j8q410f 


52859.8 


2003-Aug-09 


ACS/WFC1-2K 


1100.0 


1322.68 


1047.20 


-94.8 


60 




3(5 


j9fx01c 


53619.7 


2005-Sep-07 


ACS/WFC1-2K 


1150.0 


1317.77 


1044.48 


-95.0 


55 




37 


j9fx021 


53628.8 


2005-Sep-16 


ACS/WFC1-2K 


1150.0 


1317.91 


1044.13 


-94.8 


60 




38 


j9fx03u 


53638.7 


2005-Scp-26 


ACS/WFC1-2K 


1150.0 


1318.42 


1043.14 


-94.6 


56 




39 


j9fx04z 


53645.7 


2005-Oct-03 


ACS/WFC1-2K 


1150.0 


1318.45 


1042.15 


-94.4 


60 




10 


j9fx051 


53655.7 


2005-Oct-13 


ACS/WFC1-2K 


975.0 


1318.91 


1041.79 


-94.2 


60 




11 


j9fx06s 


53665.7 


2005-Oct-23 


ACS/WFC1-2K 


1150.0 


1319.41 


1039.97 


-93.9 


58 




42 


j9fx07x 


53673.7 


2005-Oct-31 


ACS/WFC1-2K 


1150.0 


1318.87 


1038.25 


-93.6 


52 




43 


j9fx08f 


53682.3 


2005-Nov-08 


ACS/WFC1-2K 


1150.0 


1318.60 


1036.49 


-93.2 


55 




4-1 


j9fx09j 


53690.7 


2005-Nov-17 


ACS/WFC1-2K 


1150.0 


1319.01 


1033.36 


-92.6 


56 




45 


j9fxl0u 


53699.7 


2005-Nov-26 


ACS/WFC1-2K 


1150.0 


1318.17 


1027.05 


-91.6 


53 




46 


j9fxllc 


53709.5 


2005-Dec-05 


ACS/WFC1-2K 


1150.0 


1281.89 


871.96 


-62.2 


41 




47 


j9fxl2h 


53718.6 


2005-Dcc-15 


ACS/WFC1-2K 


1150.0 


1269.60 


851.69 


-57.2 


48 





Note. — Each pair consists of two identical observations taken for cosmic-ray rejection (CRSPLIT— 2). All of the WFPC2 obser- 
vations were taken with the F547M filter, while the ACS observations were taken with the F550M filter. We processed each of the 
exposures separately. 

a Root name of the dataset in the STScI archive.* 5 Raw pixel position of the Crab pulsar. c Number of reference stars that we used 
on each image, excluding the Crab pulsar. Group number assig ned by INR061 



WFPC2 exposure-time calculat or (ETC) and sy n phot 6 , 
with the input spectrum from ISollerman et al.l (|2000f) . 
we would expect about 600 counts s -1 , or 600,000 total 
counts from the pulsar in a typical 1000-s exposure. Even 
with a gain of 15e~/ADU (used in some of the obser- 
vations, where ADU is analog-digital units), this is still 
40,000 ADU over just a few pixels, and WFPC2 saturates 
near 3500 ADU (depending on the gain). The pulsar is 
saturated in all of the F547M data by up to a factor 
of 10 in both PC an d WF observations, but this is not 
mentioned by CM99, although they do check for satura- 



tion among the reference stars. Saturation can severely 
degrade WFPC2 data since the point-spread function 
(PSF) is undersampled by the detector and most of the 
flux is concentrated in a small number of pixels. Measur- 
ing positions is particularly difficult for saturated stars, 
since most of the positional information in normal cases 
comes from the central pixels where the PSF is changing 
most rapidly. When these central pixels are saturated, 
one is forced to fit using the gently sloping wings of the 
PSF, and they provide a very weak handle on the posi- 
tion. In the next section, we discuss in detail how much 



the sat uration degrades our astrometric accuracy. In ad- 
1 See |http : //«w« . stsci . edu/resources/sof tware.hardware/stsdas /sijjifljfofi. there is no dithering between the exposures taken 
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at the same epoch; such ditherin g can help overcome th e 
undersampling of WFPC 2 fsee [Anderson fe King||2006T ). 

The a nalysis done by INR06I improved upon that of 
ICM991 [NR06 did not resample the data (resampling 
can degrade the astrometry and introduce numerical ar- 
tifacts) but treated each position measurement individ- 
ually. Additionally, they used 15 reference stars instead 
of four, used an improved distortion solution, fit for the 
proper motions of reference stars and for the orienta- 
tions and scales of the exposures, and attempted to ac- 
count for the saturation of the Crab pulsar. Finally, 
they used many more exposures. H owever , this approach 
was still not ideal, as mentioned bv INR06I themselves, as 
they used Gaussian fits instead of effective PSF (ePSF; 
lAnderson fc King] l2000f ) fits for the position measure- 
ments. Overall, including estimates of the uncertainties 
due to residual distortion error. [NR.06I find uncertainties 
of ±0.8 masyr -1 in each coordinate of the proper motion. 

We h ave attemp ted to improve on the analyses of both 
CM9l and |NR06. The most significant improvement 
comes from using many more observations: in addition 
to the large number of WFPC2 observations used by 
NR06, we used a sizable number of exposures with the 
ACS/Wide Field Camera ( ACS/W FC), a small number 
of which were discussed by |NR06| but not incorporated 
into their final analysis. We also use proper ePSF mea- 
surements and the latest distortion solutions. 

As we show, our analysis yields measurement uncer- 
tainties on the proper motion of ±0.4masyr _1 . At this 
level, we must consider in detail the reference frame of 
the measurement that we make and the corrections nec- 
essary to transform our measurement into the reference 
frame of the Crab nebula. We do this in § 12. 5[ and find 
that the reference frame uncertainties dominate the mea- 
surement uncertainties by a wide margin. 

From the WFPC2 observations, we selected only those 
taken with the F547M filter: there were too few observa- 
tions with the other filters (< 6 in a given filter) to allow 
us to properly characterize the data. We also restricted 
our data to observations where the Crab pulsar was ei- 
ther on the Planetary Camera chip (PC) or Wide Field 
chip #3 (WF3) and only analyzed the chip that the pul- 
sar was on, as for the other data sets either the pulsar 
was too close to the central reflecting pyramid 7 for reli- 
able astrometry (pixel values x or y < 100, generally) or 
there were again too few observations for proper charac- 
terization. We included only the ACS observations taken 
with the F550M filter (similar to the F547M filter), as the 
data with other filters were either too sparse to charac- 
terize or had no reliable distortion solutions/point spread 
functions (e.g., data taken through polarizers). These ob- 
servations used the WFC1-2K mode, where only one of the 
two WFC detectors is active, and only a 2048 x 2048 pixel 
sub-region of that detector is read out (of the complete 
4096 x 2048 pixel detector), thus giving a field-of-view 
that is one quarter the area of the complete ACS/WFC. 
Our final set of observations is listed in Table [I] and con- 
sists of 47 pairs of exposures, where each pair consists of 
two exposures at the same position taken for cosmic-ray 
rejection (i.e. CRSPLIT= 2). 

We took the pipeline processed images from the HST 
archive, leaving them at the flatfielded stage but not 



applying any drizzling (|Koekemoer "et~aLl [III). We 
identified cosmic rays from the CRSPLIT pairs, using 
the task driz_cr from th e STSDAS dither package 
(|Fruchter fc Mutchlerl [l99l) for the WFPC2 data and 
using the pipeline-produced data-quality extensions for 
the ACS data. For each individual exposure (we did not 
combine CRSPLIT pairs or different WFPC2 detectors), 
we performed ePSF astrometry, using the distortion so- 
lution s and ePSFs for WFPC2 from Anderson fc King 
(2003) and for the ACS/WFC from lAnderson fc Kins 
(2006). We note that the pulsar was substantially sat- 
urated on all of the exposures, both ACS and WFPC2, 
as were a small number of other stars. For the WFPC2 
data, where the undersampling is particularly bad, we fit 
the saturated stars with a larger ePSF (6-pixel radius) 
that extends further into the wings, which significantly 
improved the reliability of the measurements in tests that 
we did. For the ACS data the effects of saturation were 
not as bad, as the instrument has a larger dynamic range 
and the better spatial sampling means that more unsatu- 
rated pixels are available in th e wings of the PSF. There- 
fore we applied the standard lAnderson fc King] (|2006l ) 
ePSF technique to the ACS observations of saturated 
sources. To avoid measurements that were contaminated 
by cosmic rays, we rejected individual star positions of 
non-saturated sources if there was even one pixel con- 
taminated by a cosmic ray (identified by the algorithms 
above) within the central 5x5 pixel box used for the 
astrometry; however, we did not reject measurements of 
saturated sources (including the pulsar), since the cosmic 
ray identification routines could not distinguish between 
real cosmic rays and the effects of saturation. We also 
rejected all WFPC2 measurements with x or y < 100 to 
avoid the effects of the central reflecting pyramid. 

We assembled all of the position measurements, start- 
ing with the ACS observations. These had a 100" x 100" 
field of view and we were able to identify up to 73 stars 
besides the pulsar on those images, but we rejected 7 of 
them as they were too close to the edges. This left us 
with 66 unique stars, of which we detected up to 60 on 
any individual image. We then identified those stars on 
the WF3 and PC exposures and found that there were 
an additional 8 stars that we could identify that were not 
on the ACS images, so we have a total of 74 st ars tha t we 
used. Note that we did not use star #0 from NR06| as it 
was too saturated to have reliable measurements, but we 
have a sufficiently large number of other stars that our 
analys is is still robust. Also, the preferred solution from 
NR06] used groups 1, 3, and 6 in their numbering; their 
group 1 corresponds to our pairs 1 and 2, their group 3 
corresponds to our pairs 11-15, and th eir grou p 6 corre- 
spon ds to o ur pairs 24-33 (see Tab.[T]); 1CM99I used data 
from iNROa 's group 3 as well as two observati ons wh ere 
the pulsar was on the WF2 detector, but like INR06I we 
chose not to analyze those observations since there were 
too few to understand the uncertainties (see below). 

2.1. Measurement Uncertainty Estimation 

To properly combine all of the position measurements 
in a statistically meaningful analysis, we need estimates 
of the individual astrometric uncertainties. We took 
advantage of the CRSPLIT pairs, between which there 



should be at most a very small shift/transformation 
' See |htt P : //www . stsci . edu/instruments/wf P c2/Wf P c2_dhb/wf pc2. ch^fttn^ telescope jitter (variations in telescope pointing) 
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and breathing (variations in the detector scale d ue to 
thermal fluctuations; e.g.. lAnderson fc Kind l2006). For 
each instrument separately (ACS/WFC, WFPC2/PC, 
and WFPC2/WF3) we compared the position of each 
star with that in the other CRSPLIT image, measuring 
position differences Ax StP and Ay Si p, where s is an index 
that runs over the number of stars and p is an index that 
runs over the number of exposure pairs N p . We then 
determined for each star the variance of those position 
differences: 

N p 

°i, s = N _ 1 Ax ^p w 

" p— 1 

and the same for ov s . Without dithering, there could be 
additional uncertainties due to pixel-phase errors (errors 
from stars landing at d ifferent positions within a pixel; 
lAnderson fc Kina T2000) or from uncorrected distortion 
that we do not see from this CRSPLIT analysis, but as we 
see later our estimated uncertainties were largely suffi- 
cient. 

For the ACS/WFC data, we started with the relation 
of positional uncertainty as a function of instrumental 
magnitude (—2.5 logm Counts i n a 5 x 5 pixel box) for the 
WFC from lAnderson fe King! (pfM Fig. 13). Since the 
brightest non-saturated stars will have more un c ertain 
measurements than those in lAnderson fc King] (2006) 
and we have not derived an updated ePSF or distortion 
solution, we had an artificial minimum at 0.015 pixels for 
bright stars. To account for saturation, which occurs at 
m-inst < — 14 for the WFC, we increased the uncertainty 
to 0.075 pixels. This trend gives a reasonably good fit 
to the measured standard deviations (divided by y/2) , as 
shown in Figure [TJ and (§ I2.3[) also works well for the 
final analysis. We used several saturated stars in addi- 
tion to the Crab pulsar: they contribute very little to the 
actual fit, but by examining their residuals (e.g., Fig. [I} 
we gain a check on how well we can expect the pulsar 
data to fit. On the faint end, we included stars down to 
a WFC instrumental magnitude of —7.5: we could have 
chosen a brighter limit with smaller uncertainties, but as 
shown in the top panel of Figure [1] the number of stars 
is increasing and this increases the reliability of the fit. 
This is especially true for the WFPC2 data, where there 
are fewer reference stars overall. 

For the WFPC2 data, we started with the tr end found 
in iKaplan. van Kerkwiik. fc Anderson] |2002), and we 
used the same trend for uncertainty in pixels as a func- 
tion of magnitude for both the PC and WF3. We 
then followed the procedure outlined above, although we 
found that we had to multiply the trend for the WF3 data 
by a factor 8 of 1.5. For the saturated stars, which had 
m inst < — 10, we increased the uncertainty to 0.1 pixel 
for the PC and 0.15 pixel for WF3. 

For both the ACS and WFPC2 data, we used the fits 
for the uncertainties as a function of magnitude rather 
than the actual uncertainty for each star as the measured 
uncertainty is estimated from only a few measurements 
and is therefore noi sy, while the fit is more predictable. 
IKaplan et al.l (|2007l ) experimented in detail with differ- 
ent types of uncertainties and found that they largely did 

8 Since the trend was derived for the PC it is not surprising that 
its absolute scale should be different for the WF chips, but the 
shape seems consistent. 
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Fig. 1. — Bottom panel: Position uncertainty for each star as 
a function of instrumental magnitude (— 2.51og 10 Counts in a 5 X 
5 pixel box) for the ACS/WFC data. The crosses are the x position 
residuals, and open circles are the y residuals, and the filled squares 
are the standard deviations of x and y combined, all divided by 
V2 so as to be appropriate for a single obs ervation. The solid 
line is the uncertainty model derived in § 12.11 The Crab pulsar is 
the object at instrumental magnitude of —15.8 and our primary 
reference star (#4) is the object at —14.1, both with a red stars. 
Top panel: number of stars used as a function of m,p550M, where 
we have taken mpssoM = *"inst + 32.6. The abscissae of the two 
panels are aligned. 

not affect the final results. 

2.2. Near-IR Color-Magnitude Diagram 

The interpretation of our results depends critically on 
the distances of the field stars to which we reference 
the Crab pulsar's proper motion. We have thus de- 
rived a color-magnitude diagram of stars near the pulsar, 
using near-IR photometric obse rvations with the W ide 
Field Infrared Camera (WIRC; IWilson et all l2003f) on 
the Palomar 200-inch telescope. The observations were 
on 2003 November 19, and we exposed for 15 x 20 s in the 
J and K s filters. The seeing was not very good, and aver- 
aged 1'.'5. For the reduction, we subtracted dark frames, 
then produced a sky frame for subtraction by taking a 
sliding box-car window of 4 exposures on either side of 
a reference exposure. We then added the exposures to- 
gether, identified all the stars, and produced masks for 
the stars that were used to improve the sky frames in 
a second round of sky subtraction. We referenced the 
astrometry and photometry to the Tw o Micron All Sky 
Survey (2MASS: ISkrutskie et al.ll200l . using > 200 well- 
detected stars that were not knots of nebulosity . 

We used sextractor (Bc rtin fc Arn outs 199y) to pro- 
duce a source list for the -ftT s -band image, and then ran 
sextractor again on the J-band image using the K s - 
band source list as a set of starting positions. Finally, 
we produced the color-magnitude diagram shown in Fig- 
ure n 

Given the depth of the images and the poor seeing, 
we could not measure near-IR magnitudes for most of 
the stars that we used for astrometry. These stars were 
mostly within the Crab nebula and the high background 
limited the ground-based image. Therefore our color- 
magnitude diagram predominantly includes stars from 
the full 8' x 8' field outside the Crab nebula. We do 
not expect that the different astrometric and photomet- 
ric samples will be biased relative to one another, except 
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Fig. 2. — Near-IR color-magnitude diagram of the Crab pul- 
sar field. We plot K s magnitude vs. J — K s color for the stars 
detected on our 8' X 8' WIRC image. The uncertainties are domi- 
nated by uncertaintie s in our pho tometric zeropoints of si 0.03 mag. 
The tracks are from Cox (2000): the solid track is the main se- 
quence (luminosity class V) for Ay = 1.0 and distance lkpc, the 
dashed track is the main sequence for Ay = 1.0 and d = 4kpc, 
and the dot-dashed track is the giant branch (luminosity class III) 
for Ay = 1.0 and d = lOkpc; some stellar types are labeled. We 
also show a reddening vector for Ay = 2.0. The diagonal lines 
are approximate lines of constant V (or F550M) magnitude, as- 
suming V — K B Si 4.3(J — K B ). Note that V magnitude is very 
close to both mp547M an d "Ifssom- The stars are the 10 objects 
used in our HST astrometry that we could detect in the WIRC 
images, the Cra b pulsar is the green diamond (labeled), and star 
#4 from NR06 (which is our proper motion reference star) is the 
green square. 

that the limiting magnitude in the near-IR is brighter: 
the photometric stars are those that are either very 
bright and/or are outside the Crab nebula, while the 
astrometric stars include those inside the nebula for a 
range of brightnesses. We estimated a rough conversion 
between our near-IR photometry and our HST photom- 
etry using synphot for Ay — 1.0 mag, where we find 
V - K s « 4.3( J - K a ) and V « m F 550M ~ m F5i7M 
(since the STMAG system is based on the Vega system 
at V, no zeropoint offset is necessary). A rough photo- 
metric calibration for the HST data can be done with 
the instrumental magnitudes that we measure, such that 
TOF550M ~ m-inst + 32.6, where we neglect any aperture 
corrections. As can be seen from Figure [1] the major- 
ity of the stars are at tofssom ^ 21, which is fainter 
than most stars in Figure [2] but does not imply drasti- 
cally different distances. The WFC saturation limit of 
wiinst = — 14 translates to TOF550M ~ 19, and indeed the 
stars above this line in Figure [2] are saturated. 

From this we see that the majority of the stars in this 
field are consistent with main-sequence stars at distances 
of 1-4 kpc; a few may be more distant giants, and some of 
the astrometric reference stars may be among these, but 
they could also be slightly more reddened main-sequence 
stars (note that we are not accounting for the effects 
of metallicity). There should not be many stars that are 
much closer, as they would have to be rather redder than 
the more distant stars that we measure. 

2.3. Fitting for the Astrometric Parameters 

Once we have our data set with uncertainties as derived 
in § 12.11 we are in a position to fit for the positions and 
proper motions of each star, as well as the plate-scales, 



orientations, and c entral position of each exposure (see 
iKaplan et al.ll2007l for a detailed description of the fitting- 
procedure). For the fitting, we need to set the absolute 
plate-scale and orientation, which we did by assuming the 
plate-scale and orientation from the first exposure from 
pair 35 are the nominal values of 50masyr -1 and the 
rotation from the image header; this exposure is no more 
likely than any other to have the correct plate-scale or 
orientation, but that will just lead to absolute uncertain- 
ties on the proper motion of 0.1% or less 9 . We then also 
assumed that the orientation of the first exposure from 
pair 36 is the header value, since without two exposures 
with known orientations the fitting procedure could lead 
to a net rotation of the data with time that is compen- 
sated by a bulk proper motion. The fitting process might 
also compensate for a secularly increasing offset between 
exposures by introducing a fictitious bulk proper motion 
for the ensemble of stars. In order to prevent this, we ini- 
tially assigned a star to have zero proper motion. This 
does not actually define the reference frame, since we can 
arbitrarily shift the proper motio ns of a ll of the stars. For 
this fixed star we chose #4 from INR061 : this star has the 
advantage that it is close to the Crab pulsar, so it is on 
almost all exposures (73 of 94, after rejecting individual 
measurements for cosmic rays as described above) and 
is bright but not saturated (see Fig. [3]). From § 12.21 and 
Figure [2] star #4 appears to be roughly at ~ 4 kpc. 

Once we have measured source positions and estimated 
uncertainties at each epoch, we directly use the different 
observations to measure a proper motion. We fit simul- 
taneously for the positions and proper motions of each 
star (with the proper motion of #4 fixed to zero) and 
the transformation parameters for each exposure (with 
the rotations and plate-scales fixed as discussed above 
for the first exposures in pairs 35 and 36), including all 
of the WFPC2 and ACS data in the fit. Each exposure 
had a 6-parameter t rans formation, such as that u sed by 
IKaplan et all (|2007l ) and I Anderson fc Kind (|2006f) . Such 
a transformation is able to deal implicitly with linear 
variations in th e distortion cau sed by breathing or sys- 
tematic effects (| Anderson! |2007( ) . For the ACS data the 
fits give scale uncertainties of « 0.0005% and position 
angle uncertainties of 0?0003, and shift uncertainties of 
< 0.01 pixel. For the WFPC2 data the results are some- 
what worse, largely due to the smaller number of stars: 
WF3 has scale uncertainties of ss 0.02% and position an- 
gle uncertainties of 0?005, and shift uncertainties of 0.02- 
0.1 pixel (depending on the number of stars included), 
while the PC has scale uncertainties of ss 0.07%, rotation 
uncertainties of 0?05, and shift uncertainties of 0.1 pixel. 

Initially we achieved a reasonable fit, with \ 2 = 5935.9 
for 3808 degrees of freedom (dof; we had 2322 obser- 
vations of both x and y for 4644 data points, and 836 
free parameters), or Xred = 1-559. This x 2 results from 
comparing the computed positions of every star at every 
epoch (based on the fitted reference positions, proper 
motions, and frame transforma tions) to the measured 
positions; see Ka plan et alJ (|2007f ). Eqn. A6. However, 
there were anomalously large contributions to the total 

9 We have verified this by fitting our final reference positions 
to positions derived from the WIRC data, which are referenced to 
2MASS, which is tied to the International Coordinate Reference 
System. We fi nd r otations of < 0.1° and scal e changes of < 0.1%. 
Also sec § 12.41 and van dcr Marcl et al. ( 2003). 
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Fig. 3. — Proper motion of the Crab pulsar, shown on a drizzled lIKoekemoer et ail [20021 ACS/WFC image (pair 35). The vector 
indicated by the solid line shows the proper motion of the pulsar, corrected according to § [23] for the effects of differential Galactic rotation 
and solar motion, and showing the position of the pulsar fOOO yrs from now. The inner-most circle is the l-a statistical uncertainty; the 
middle circle is the l-cr combined (measurement error plus correction for DGR and LSR plus inclusion of uncertainty in the progenitor's 
peculiar motion relative to its lo cal stan dard of rest) uncertainty; and the outer circle the is 3-c combined uncertain ty. We also label our 
primary reference star (#4 from NROrf) and the nebula's symmetry axis (at position angle 304° north through east; Ng & Romani 2004) 
which we take to be the projected rotation axis of the pulsar. The magenta dot-dashed line indi cates ou r proper motion projected back to 
1054 CE and the green circle indicates the divergent point (with ±1-<t uncertainties) found by WM77 for the Crab's filaments. North is 
up, and east to the left. 



\ 2 from a few deviant data points. Beyond the cosmic- 
ray rejections discussed above, we rejected an additional 
8 measurements that deviated by more than 10a from 
the best-fit model, and reduced \ 2 to 4473.7 for 3792 dof 
(Xred = 1-180; the proper motion of the pulsar changed 
by < In after the rejections). The rejected measure- 
ments were distributed among the stars and exposures, 
and likely represented statistical fluctuations or unde- 
tected cosmic-rays. After the additional rejections, the 
fit looked good overall, with no individual star or expo- 
sure dominating the fit. Note that the pulsar is satu- 
rated, and so its astrometric position uncertainty is sig- 
nificantly higher than most of the reference stars. As 
such, it does not dominate the overall fit. We tested us- 
ing other reference stars and exposures (both ACS and 
WFPC2 observations), and the results did not depend on 
those choices except for a net shift in the proper motion, 
but we correct for this below. 

As discussed above, all of the proper m otions that we 
fit for are relative to that of star #4 from lNR06l . but we 
of course do not know what the proper motion of star #4 
is, and it does not make a useful reference frame: with a 
transverse velocity dispersion of ~ 20kms _1 and a dis- 



tance of a few kpc, the proper motions of random stars 
are ~ lmasyr -1 . We must therefore try to determine a 
reference frame for our measurements in which the pro- 
jected motion of the pulsar can be compared sensibly 
with the projected orientation of the nebular symmetry 
axis. We do so in two steps. First, we determine the av- 
erage proper motion of the ensemble of reference stars. 
Next, we consider how those stars are moving relative to 
the Sun. 

In Figure 2] we plot all of the proper motions that we 
measure. The Crab pulsar clearly has a much larger and 
more significant proper motion than the other objects, al- 
though we find a number of other stars with > 3er detec- 
tions of proper motion. In that figure we have determined 
the mean proper motion of all of the stars excluding the 
Crab pulsar (iteratively rejecting outliers) and shifted the 
proper motions to have zero mean. This shift has a mag- 
nitude of (Ajtia, Afig) = (-1.8±0.2,-0.0±0.2)masyr -1 , 
moving star #4 to its position away from the origin. The 
circle in Figure 0] shows the standard deviation of the 
proper motions, and it is comparable to the magnitude 
of the shift. The uncertainty in the shift is much smaller, 
though, since it is the standard deviation divided by the 
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Fig. 4. — Proper motions for each star, shifted so that the net 
proper motion is zero (represented by the dotted lines). The Crab 
(green diamond) and star #4 from NR06 (green square) arc identi- 
fied. The red ellipse shows the standard deviation of the reference 
stars (after rejecting outliers), and the magnitude of the shift is 
shown by the resulting proper motion of star #4 (whose proper 
motion was fixed to zero during the fit). 

square root of the number of stars used (here 49), and 
the shift in Right Ascension, at least, is statistically sig- 
nificant. 

To see how our results depended on the choices of refer- 
ence star/epochs, and on the datasets we used, we made a 
number of different fits. Using the whole ACS+WFPC2 
dataset we iterated among a variety of reference stars, 
choosing ones that were relatively bright but not satu- 
rated and were close to the pulsar. We also used different 
ACS reference exposures and found that the proper mo- 
tion values (corrected to have zero net stellar proper mo- 
tion) changed by at most 0.6 mas yr -1 , consistent with 
our uncertainties given below in Equation [2j 

We then tried to choose different data sets, restrict- 
ing ourselves to only the ACS data, only th e WFPC2 
data, only the WFPC2 data used by [NR06, and some 
other combinations. In particular, we fit using none of 
the stars with instrumental magnitudes (from the ACS 
data) fainter than —9.5: this excludes the portion of Fig- 
ure[T]where the trend of uncertainty vs. magnitude climbs 
upwards. We also fit excluding all of the saturated stars 
from the ACS data, with the exception of the Crab pulsar 
(of course): this fit required us to remove some of the PC 
epochs, as without the saturated stars there were too few 
reference stars for a constrained fit. Such fits proceeded 
in exactly the same manner as the fit described above, 
with the only difference being that we used different ob- 
servations as the initial reference; no other special ma- 
nipulation was required for these fits. We found that the 
corrected proper motions from all of these data sets were 
consistent with each other, as listed in Tableland shown 
in Figure El but were not necessarily consistent with the 
value s from the literature. Our result using just the data 
from iNROa was consistent with their result, although it 
had significantly higher uncertainties. This is because we 
took an uncertainty of 0.15 pix for the WF3 obser vations 
of the pulsar, which is 15 mas, while Figure 3 of[NR06 
shows their individual data points as having uncertain- 
ties of 5-8 mas. In contrast, our measurements using a 
longer time baseline and with the highe r-precision ACS 
data were not consistent with the INR06I values. In all of 
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Fig. 5. — A plot of all of the measured proper motions for 
the Crab pulsar from this work and the literature, as given in 
Table [2] these values are the direct fitted results, before correction 
for Solar motion or Galactic rotation. The unfilled triangles (at 
various orientations) are those values from the literature, while the 
other symbols are our values using various subsets of the data [the 
"WF3(NR06)" value is our measurement using just the WF3 data 
used bv lNR06|| . The best measurement is the "ACS+PC+WF3" 
value, but all of our measurements are consistent with each other. 

these fits the resulting x 2 values were close to 1.0, indi- 
cating that our uncertainty estimates (§ 12. ip — which we 
derived just by comparing pairs of observations — were 
reasonable for the dataset as a whole, and gave good val- 
ues for the different instruments. The fits that excluded 
the faint or the saturated reference stars could be used as 
our "default" fits, but we chose to retain as many stars 
as possible. With the exception of the PC epochs, where 
there are few stars, the saturated stars largely ride along 
with the fit and contribute very little, but they allow us 
to examine the goodness-of-fit for saturated objects be- 
sides the Crab pulsar, and indeed we find that the data 
fit them reasonably well (reduced x 2 of 0.85 to 1.48). 
The fainter stars contribute more to the fit, but are not 
dominant, and they allow us to examine the quality of 
the fit for a larger number of objects (again, it is good). 

Using all of the data and our standard choices for refer- 
ence exposures/star, the proper motion for just the Crab 
pulsar was a good fit, with % 2 ed = 217.3/184 = 1.18, al- 
though this x 2 lS n °t formally correct as the parameters 
for the exposures are not properly counted. We show the 
astrometry in Figure El From this one can see that some 
pairs of exposures seem to deviate systematically from 
the overall trend, such as some of the WFPC2 data from 
near MJD 51,800. These are likely related to the changes 
in position angle of the observations (Tab. Q}, although 
whether it is an intrinsic effect of the position angle (i.e. 
uncorrect ed distortion, perhaps related to charge transfer 
efficiency; iKozhurina-Platais. Goudfrooij. fc Puzial I2007T ) 
or something to do with the changing set of reference 
stars, we cannot determine. Overall, though, the devia- 
tions are not greatly significant, and the proper motion is 
confirmed by our analyses of various subsets of the data. 
The corrected proper motion is: 

fi a = -12.0 ± 0.4 mas yr -1 

H5= +4.1 ± 0.4 mas yr" 1 , (2) 

where the uncertainty is a combination of the uncertainty 
on the derived proper motion and the uncertainty in the 
shift to the reference frame. This proper motion is in 
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TABLE 2 

Measured Proper Motions for the Crab Pulsar 



Reference 


Data Source 


Data Subsets 


Ma AM 








(mas yr" 1 ) 



WM77 Plates 

CM99 Limited HST WFPC2 

NR06 HST WFPC2 

NR06 HST WFPC2 

This work HST WFPC2 

This work HST WFPC2 

This work HST ACS 

This work HST ACS+WFPC2 

This work HST ACS+WFPC2 

This work HST ACS+WFPC2 

This work HST ACS+WFPC2 

This work HST ACS+WFPC2 



NR06 groups 1, 
NR06 groups 3, 
All i'< u I .; 

NR06 groups 1, 



3, 6 
6 

3, 6 



All PC 
All WF3 

No faint stars a 

No saturated stars b 



-13 ±2 

-17 ±3 

-15 ±0.8 

-10.9±2.2 

-15.6±1.8 

-13.0±2.4 

-9.8±1.7 

-12.5±0.6 

-12.1±0.6 

-12.0±0.5 

-12.3±0.4 

-12.0±0.4 



+7 ±3 

+7 ±3 

+1.3±0.8 

+1.0±2.0 

+3.7±1.7 

+4.6±2.3 

+4.2±1.8 

+4.1±0.6 

+4.5±0.6 

+3.9±0.5 

+4.2±0.4 

+4.1±0.4 



Note. — All proper motions arc the values from fitting before correction for Solar motion 
or Galactic rotation. Proper motions from this wor k and|WM771 are explicitly in a reference 
frame defined by the background stars in this field; CM99 explicitly assumes that the stars 
have zero proper motion. Also see Figure [5] 

a We rejected the 27 stars with mi nst F550M > —9.5 (the upward part of the trend in 
Fig.^.k We rejected the four saturated reference stars. This also required removing several 
of the PC epochs since they had too few stars for proper solutions. 
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Fig. 6. — Proper motion of the Crab pulsar. In the top two panels 
we show the position in Declination (top) and Right Ascension 
(middle), plotted against MJD; Act = and AS = correspond to 
our reference epoch of pair 35 (MJD 52859.8). The observations 
with ACS/WFC are the blue stars, those with WFPC2/PC are 
red circles, and those with WFPC2/WF3 are green diamonds (as 
labeled). The magenta lines show our best-fit proper motion along 
with ±lcr uncertainties. The bottom panel shows the position in 
Right Ascension versus that in Declination: black lines connect 
each observation to the point on the best-fit proper motion line 
at the time of that observation. We als o plot t he best-fit proper 
mot ions alo ng with ±1<t uncertainties for lCM9"9l (dot-dashed lines) 
and INR06I (dashed lines); all proper motions are the final fitted 
results before correction for Solar motion or Galactic rotation. 



a reference frame denned by the average motion of the 
background stars in the field (this is th e same proce- 
dure used for the initial measurement of IWM77I . which 
is their Equation 2, but the different choices of references 
stars means that the reference frames will not be exactly 
the same). However, as we discuss below, this reference 
frame is not the appropriate one for considering the de- 
gree of alignment between the pulsar motion and the 
projected symmetry axis of the surrounding nebula. 



2.4. Absolute Position 

We can get a position for the Crab pulsar tied to 
the International Celestial Reference System (ICRS; 
useful for analysis of timing data, for example) di- 
rectly from 2MASS, which is tied to the ICRS 10 
to better than O'.'l. The pulsar is listed as source 
2MASS J05343194+2200521, at position (J2000): 

a = 05 h 34 m 31!94 

5 = +22°00'52"l, (3) 

with quoted uncertainties of ±0'/06 on each coordinate. 
This position has been precessed to equinox J2000, but 
it was actually measured on 1997 October 18. The 2.2 
years between those dates implies a shift due to the mea- 
sured proper motion of 27 mas, which is smaller than 
the measurement uncertainty. To verify the tie of the 
2MASS data to the ICRS, we compared the position of 
32 stars located within 5' of the pulsar that were detected 
in both 2MASS and the Second U. S. Naval Observa- 
tory CCD Astrograph Catalog (UCAC2: IZacharias et al.l 
I2004D . which is tied directly to the ICRS. We find a small 
residual shift (« 32 mas) between 2MASS and UCAC2 
(the UCAC2 proper motions for these stars are not sta- 
tistically significant, so we did not include them), but 
this is again smaller than the measurement error re- 
ported by 2MASS. The 2MA SS position is con sistent 
with the radio imaging position (|Han fc T ian 1999) listed 
in SIMBAD from t he NRAO VLA Sky Survey (NVSS; 
iCondon et al.lfi998h . but should be more accurate. A 
caveat, though, is that 2MASS could not measure the 
pulsar sep arately from the op tical/infrared knot located 
0'.'7 away (H ester et al.l Il995h . However, the pulsar is 
> 40 times brighter than the knot (at least at 5500 A), 
and so the 2MASS centroid should essentially be at the 
pulsar's position. 

2.5. Proper Motion Reference Frame 

10 See http: //spider . ipac . caltech. edu/staf f /hlm/2mass/overv/overv.htm 
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Equation [2] gives the proper motion of the Crab 
pulsar relative to an ensemble of stars. Normally, 
we would correct our measured proper motion for 
the peculiar velocity of the Sun relative to the lo- 
cal standard of rest (LSR; we take the Solar motion 
to be [U,V,W] e = [10.00,5.25,7.17] kms" 1 , where 
the uncertainties on those components are negligible 
compared to the other uncertainties discussed below; 
iDehnen fc Binnevl [19981 and for the effects of differen- 
tial Galactic rotation (PGR; we use th e Galactic poten- 
tial model of iKuiiken fc Gilmorei [l989l to determine the 
rotation curve, and then ex amine deviations from the ro- 
tation curve determined bv lBrand fc Blitz 1993|), so that 
the proper motion reflects a velocity relative to the ob- 
ject's local standard of rest. However, we should not 
apply these corrections blindly for the Crab pulsar. The 
reason is that the astrometric reference stars that we are 
using are at distances of 1-4 kpc (and possibly a bit fur- 
ther), comparable with the nominal 2 kpc distance of the 
Crab pulsar 12. 2p . To first order, then, we do not need 
to correct for DGR or LSR motion. 

To see ho w valid t his is, we refine our distance estimate 
by following|WMlI They found a mean statistical paral- 
lax of 0.5 mas for the reference stars from kinematic con- 
straints, roughly consistent with our estimate of 1-4 kpc 
distanc es for the WIRC stars. Indeed, the stars used by 
WM77 were spread over several arcminutes and there- 
fore are more compara ble to the near-IR sample than the 
HST sample. IWM77I formed their estimate by examin- 
ing the dispersion of the reference stars about their mean 
proper motion and then relating that to the expected ve- 
locity dispersion. We find a proper motion dispersion of 
o-p = 1.8 masyr . Comparing this to the expected 1- 
D velocity dispersion of S ^kms" 1 (appropriate for la te- 
type main-sequence stars; iBinnev fe Merrifieldlll998l p. 
632), we find a mean distance of 3-4 kpc, consistent with 
our photometric estimates and with the likelihood that 
some of the HST stars are fainter and more distant than 
the near-IR stars. 

Of course, the stars are not all at a single distance, 
but are distributed over a range, and their contribution 
to the ensemble's proper motion depends on their bright- 
nesses and on how many observations we have for them. 
The range of distances will influence the ensemble's mean 
proper motion: going from 1 kpc to 4 kpc increases the 
correction due to DGR by 0.2 masyr" 1 , while the correc- 
tion for the LSR decreases in magnitude by 1 masyr" 1 . 
Since the Crab pulsar is toward the anti-center the LSR 
correction is more significant, and it depends inversely 
on distance. We take the pulsar to be at 2 kpc, which 
gives corrections of 

(A/it, A Mb ) DGR (2 kpc) = (+0.68, +0.04) masyr" 1 

(Am, A/i b ) LSR (2kpc) = (+0.47, -0.64) masyr' 1 . (4) 

For the average stellar frame at 4 kpc, the corrections 
are: 

(Afi e , A^) DGR (4kpc) = (+0.78, +0.03) masyr" 1 
(Afi t , A/i b ) LSR (4kpc) = (+0.23, -0.32) masyr" 1 , (5) 
so the net corrections are 

(Afj, e , A/i fe ) D Ga(4kpc -> 2 kpc) = (-0.10, +0.01) masyr" 1 
(A/if, A^ 6 ) LSR (4 kpc -» 2 kpc) = (+0.24, -0.32) masyr" 1 . 

(6) 



In equatorial coordinates these corrections are 
(Afi a , A/x 5 ) D GR(4kpc -> 2 kpc) = (-0.05, +0.09) masyr" 1 
(A/i Q , A^) LSR (4 kpc -> 2 kpc) = (-0.14, -0.37) masyr" 1 . 

(7) 

These corrections must be subtracted from the proper 
motion we found in § 12.31 so the proper motion of the 
pulsar relative to its local standard of rest is fi a = 
— 11.8 masyr" 1 , = +4.4masyr" 1 . Of course, these 
distances are uncertain. Assuming a 0.5 kpc uncertainty 
in the distance to the pulsar, and considering the range 
of 3 kpc to 5 kpc for the reference stars, we found that 
the correction overall varied by ±0.1 masyr" 1 from the 
nominal values in Equation [71 so compared to the other 
sources of uncertainty this is a minor effect. We note 
that the uncertainties here exceed formal 1-er confidence 
intervals, since the distance intervals represent the full 
range of plausible distances. We also note that we have 
chosen a particular formulation of the Galactic rota- 
tion c urve through the potential of IKuiiken fc Gilmorei 
(1989): other formulations (flat rotation curves, use of 
Oort constants, etc.) give slightly different results for 
the DGR corrections. In particular, without a velocity 
that varies as a function of distance, the correction from 
4 kpc to 2 kpc is ident ically 0. Using other choices (e.g., 
the rotation curve of iBrand fe B litz 1993) changes the 
correction by a small amount on an absolute scale, typ- 
ically ±0.1 masyr" 1 (other determinations of the LSR 
corrections agree to within 10%). We include this as an 
additional uncertainty. 

However, more significant uncertainties come from our 
assumption that the rotation curve toward the outer 
Galaxy is both well-known and well-behaved, when it 
is neither. Overall, the rotation curve of the outer 
Galaxy has line-of-sig ht random variations of ±5kms _1 
(|Brand fe B litz 199 J and references therein), which ac- 
counts for the deviations of individual locations from the 
bulk velocity. The curve itself is poorly measured toward 
the anticenter, but this value should be relatively inde- 
pendent of position. So we would expect 1-D variations 
of ±5kms" 1 on top of the circular velocity, implying a 
proper motion uncertainty of 0.5 masyr" 1 at 2 kpc. The 
proper motion is 

fj, a = -11.8 ± 0.4 ± 0.5 masyr" 1 

fj, s = +4.4 ±0.4 ±0.5 masyr" 1 , (8) 

compared to the local standard of rest of the Crab pul- 
sar, where the first uncertainty is the measurement un- 
certainty (Eqn. [2]) and the second is from the reference 
frame uncertainties. 

Now, if our goal is to compare the proper motion vector 
to the projected orientation of the torus axis, we need to 
account for yet another zero point uncertainty, namely 
the unknown velocity of the Crab's progenitor. The pro- 
genitor was presumably not stationary with respect to 
the local standard of rest at that position, and any pecu- 
liar velocity should remain after the explosion. Therefore 
a final zero-point uncertainty comes from the unknown 
peculiar velocity of the Crab pulsar's progenitor. First, 
there are bulk streaming motions in th e outer Galaxy 
with line-of-sight mag nitudew 12kms" 1 (|Brand fe Blitzl 
119931) . Beyond this, early-t ype stars have 1-D velocit y 
dispersions of - lOkms" 1 flBinnev fc Merrifieldl [1991 . 
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Together these give into a 1-D velocity uncertainty of 
16 kms -1 , which, combined with the 0.5 masyr -1 uncer- 
tainty discussed above, translates into an uncertainty of 
1.7masyr _1 at 2kpc, with a range of 1.4-2.3 masyr -1 
for a distance range of 1.5-2.5 kpc (i.e. an uncertainty 
on the uncertainty). We assign a conservative value of 
2.0 masyr -1 (considering the uncertainty on the distance 
as well as the reference frame velocities) for the sys- 
tematic uncertainty due to combination of the reference 
frame effects, but also note that the progenitor's velocity 
could have been much larger (see below). We then find 
a proper motion in the reference frame of the progenitor 
star of: 

pL a = -11.8 ± 0.4 ± 2.0 masyr -1 

Hs= +4.4 ± 0.4 ±2.0 masyr -1 , (9) 

where again the two uncertainties are from the measure- 
ment and the unknown reference frame, as discussed 
above. This proper motion has a magnitude of fi — 
12.5 ± 0.4 ± 2.0 masyr -1 at an angle of 290° ± 2° ± 9° 
(east of north), or a transverse velocity of 120 kms -1 for 
a distance of 2 kpc; the velocity is quite uncertain, both 
due to the uncertain frame of the proper motion and the 
distance uncertainty. 

3. DISCUSSION AND CONCLUSIONS 

While we assumed a velocity of ±10 kms -1 for the 
Crab pulsar's progenitor, it is possible that the progen- 
itor itself had a mu ch larger space ve locity. This was 
proposed early on by [Minkowski (|1970t ). who considered 
the high proper motion of the pulsar to be a relic of the 
progenitor's velocity and that the progenitor itself was 



a runaway star (e.g., lBlaauwlll961) from the Gem OBI 
association. However. I Minkowskil ( 19701 ) later dismissed 
this hypothesis since he did n ot believe t h at th e super- 
nova was a type II explosion. iGott et ail (|1970f ) have a 
similar idea, where they consider the Crab pulsar and 
the nearby pulsa r B0525±21 (two of the fi rst pulsars 
to be discovered; iStaelin fc Reifensteinl [1968) as former 
binary companions ejected from the Gem OBI associ- 
ation; lHarrison et all |l993) attempted to measure the 
proper motion of PSR B0525+21 but did not find a sta- 
tistica lly sig nificant result. Later analyses , such as iPolsl 
(| 19941 ) and iMdzinarishvili fc Dzigvashvilil (|200lD . have 
revived the hypothesis that the pr ogenitor had a large 
(> 100 kms -1 ) space velocity, with iPoisl (|1994f ) arguing 
that the large height below the Galactic plane (200 pc for 
a distance of 2 kpc , which is larger than the scale-height 
of OB stars; [Reedl [2000l: lElias. Cabrera-Cafio. fc Alfarol 
2006) and the presumed evolutionary state of the pro- 
genitor (inferred from the elemental and velocity struc- 
ture of the Crab Nebula) suggest that the Crab pulsar 
was formed from the second explosion in a binary sys- 
tem, and that it had a large space velocity from the first 
explosion. Currently we cannot say definitively whether 
or not the progenitor had a significant space velocity and 
must therefore treat this as an overall uncertainty on our 
whole analysis. 

3.1. Location of the Explosion Center 

A number of authors have estimated the "divergent 
point" for the Crab Nebula: the point from which all of 
the filaments seem to traveling outwards, which is pre- 
sumed to be the center of the explosion. Among the 



more reliable measurements are those of lTrimbld (jl968ft . 
IWM771 . and iNugentl (|1998l who also review the situa- 
tion), which we list in Tabled Most of these determi- 
nations trace the filaments back and find best-fit dates 
for the explosion o f » 1130 CE instead of the commonly 
accepted 1054 CE (jSteph cnson fc Greenl l2002). with the 
difference caused by unmodeled acceleration of the fila- 
ments. 

With the proper motion that we derive, the explosion 
centers from the literature, and the nominal explosion 
date of 1054 CE, we can perform three tests: we have a 
time, a displacement, and a velocity, and we can use any 
two of those to estimate the third. First, we can mea- 
sure how close our proper motion comes to the various 
explosion centers for the nominal explosion date, and we 
give these values in the last column ("Ar(1054 CE)") of 
Table [31 Second, we can compute the dates of closest ap- 
proach between our projected proper motion vectors and 
the explosion centers, which serve as our own estimates 
of the explosion dates. The uncertainties on those val- 
ues are dominated by the uncertainties of the divergent 
point measurements (~ 1"), and the values are given in 
the "Divg. Date" column of Table [3j with the approach 
distances for those dates given in the "Ar m ; n " column. 
Finally, we can compute our own estimate for the ex- 
plosion position, taking our proper motion and assuming 
the date of 1054 CE, and we give this in the last row of 
Table 

In general, all of the values — the approach distances 
for 1054 CE, the explosion dates, and our inferred explo- 
sion center — are consistent at better than 1-cr. The first 
two elements are largely consistency checks: this shows 
that our proper motion is consistent with the indepen- 
dent divergent point estimates, and that our reference 
frame corrections are consistent (alt hough t hey were not 
explicitly the same). As discussed in lWM77l . the location 
of the divergent point depends on the choice of reference 
frame, so we cannot address any of the larger reference 
frame uncertainties. We note, though, that in contrast to 
our proper motion that approaches the divergent points 
with distances of < 1", the proper motion of INR06I does 
approximately a factor of 3 worse. The final element 
that we have computed, our own estimate for the explo- 
sion center, is more precise than previous estimates by 
a factor of ~2-3 in each axis. This may serve to help 
constrain future measurements of the filament motions 
and acceleration, as its independence from the filaments 
themselves should improve the reliability of the measure- 
ments. 

3.2. Spin-Kick Mis Alignment 

iNg fc Romanil (|2004l ) fit a model to the torus seen in 
the HST data, and find a best-fit torus symmetry axis of 
304.0° ± 0.1°, which is the projection of the spin axis on 
the plane of the sky. This implies a projected misalign- 
ment of 14° ± 2° ± 9°, as seen Figure [3] This p rojected 
misalignment is less than that found in INR06I . and is 
significant if one only considers the measurement uncer- 
tainty: with all of the uncertainties, the misalignment is 
consistent with a broad range of values, including zero. 

Perhaps the next best case of a pulsar wind nebula 
giving the projected spin axis is that of the Vela pulsar, 
where the proper motion (corrected for Galactic rotation 
and solar motion) is 45 ± 1.3 masyr -1 at a position an- 
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TABLE 3 

Comparison Between Our Proper Motion and Divergent Points 



Reference 


Divergent Point (J2000) a 


Divg. Date b 


Ar min c 


Ar(1054 CE) d 




a d 


(CE) 


(arcsec) 


(arcsec) 


Trimble f 19681 


05 h 34 m 32f 72 ± 0?12 +22°00'47'.'5 ± 1'.'4 


1067 ± 138 


0.5 


0.6 


WM77 


05 h 34 m 32!67 ± 0?06 +22°00'47'.'6 ± 0'.'9 


1114 ±78 


0.7 


1.0 


Nugent (1998) 


05 h 34 m 32!84 ± 0?12 +22°00'48'.'0 ± 1'.'3 


947 ± 138 


0.6 


1.5 


This work c 


05 h 34 m 32f 74 ± 0?03 +22°00'47'.'9 ± 0'.'4 


1054 







a The positions of the di vergent point according to [Trimble (1968) and Nugent ( 199g) were computed 
from the offsets given in [Nugent (1998) between the divergent point /explosion center and the star 5" 
to the north-east of the pulsar, whose position we tak e to be: a = 05 h 34 m 32?17, 8 = +22°00'56'.'0 from 
2MASS (the star is 2MASS J05343217+2200560). For lWM77l . we take the divergent point directly from 
their paper (Eqn. 17). b Date of closest approach between our proper motion vector projected backwards 
and the divergent point. The uncertainties are only measurement uncertainties — no reference frame 
uncertainties are included. Closest approach between our proper motion projected backward and the 
estimated divergent point of the filaments. d Distance between our proper motion projected backward 
and the divergent point for 1054 CE. e Not truly a divergent point, but rather the location of the pulsar 
projected back to 1054 CE. 



gle of 301° ± 2° (jDodson et al.l [2003h . The symmetry 
axis of the torus is a t a position angle of 310.6° ± 0.1° 
(|Ng fc Romanil[200l. giving the projected misalignment 
of 10° ± 2° (|Ng fe Romam1l2007f) . For this system, the 
proper motion reference frame and corrections should be 
better defined than those for the Crab: the proper mo- 
tion is measured in the radio, so the reference sources 
are at infinite distance; and the Vela pulsar is reason- 
ably close (with a distance measured through geometric 
parallax) and not located at the Galactic anti-center, so 
the effects of Gal actic rotation are mu ch better under- 
stood. However, iNg fc Romanil (|2004h still fail to in- 
clude any allowance for the unknown velocity of the pro- 
genitor: a ztlOkms" 1 velocity at a distance of 287 pc is 
±7masyr _1 (or an angular uncertainty of ~ 9°), so again 
this completely dominates the measurement uncertainty 
and makes the degree of misalignment consistent with 
zero. 

The two examples considered here, the Crab and Vela 
pulsa rs, while the two best X-ray tori (|Ng fc Romanil 
2004), may not be good cases for computing misalign- 
ment. This is because they are both moving at smaller 
velocities than the average pulsar pop ulation (120 km s -1 
and Glkms" 1 , vs. ~ 400 km s" 1 ; e.g. jHobbs et alj|2005t 
iFaucher-Giguere fc Kaspil2006l ). so the unknown velocity 
of the progenitor has a correspondingly greater contribu- 
tion. In fact, there is a bias in favor of tori being found 
around relatively slow pulsars. This is because the tran- 
sition from a "bubble" pulsar wi nd nebula (PWN) with a 
torus to a bow-shock PWN (see iGaensler fc Slane!l2006T ) 
occurs when the pulsar has traveled roughly 68% of the 
distance fro m the center of the of the supernova remnan t 
to its edge (|van der Swaluw, Downes, fc Keeganl l2004f ). 
largely independent of the pulsar's velocity. So slower 
pulsars will spend longer in the bubble/torus phase. For 
faster moving pulsars, whose proper motions and pro- 
jected spin axes are not as well determined (in general 
they are further away), the uncertainty will be closer 
to 1° and will not dominate over the me asurement un- 
certai nties. In addition, if the model of INg fc Romanil 
(2007) is correct, we would expect the faster moving pul- 
sars to be intrinsically closer to alignment. We also note 
that, for pulsars moving at > 200 kms -1 , all of the refer- 



ence frame uncertainties discussed here will lead to un- 
certainties of < 10% on the space velocity: this will typi- 
cally be less than the uncertainty on the distance. There- 
fore, in studying the magnitude of pulsar velocities (or of 
the pulsar population) the reference frame uncertainties 
will not be significant. 

However, we can also look at the situation from the 
other side. If the spin and kick axes were perfectly 
aligned in the reference frame of the progenitor's mo- 
tion, we could still (erroneously) infer a misalignment 
because of a high progenitor space velocity. In practice, 
though, it is difficult to disentangle the effects of intrinsic 
and apparent misalignments for single objects. A further 
complication comes from the fact that all alignments are 
examined only in projection: inclinations of nebulae can 
be estimated (although not directly measured), but ra- 
dial velocities of pulsars are entirely unknown, and with- 
out the third dimension any observed alignments could 
still be coincidences. 

A number of o ther pulsar wind nebulae also have sym- 
metry axes je.g jPavlov et alJl200H : iHelfand et alJl200H : 
Ng & Romani 2004) , although the lower fluxes and larger 
distances make most of these hard to observe in de- 
tail. There are other situations where alignments are 
inferred but not measured directly: for instance, using 
an offset from the center of a supernova remnant to de- 
rive a kick direction (this approach can introduce sub- 
stant ial systematic errors of its own; sec Gaensler et al. 
2006) , and deriv ing a rotation axis from fitting radio po- 
larization data dDeshpande et all 119991: lLai et all 120011 : 
iRomani fc Nsll2003l : lWang et al 1120061 : lRankinll2007h ■ For 
those syste ms with constrained proper motions and rota- 
tions axes. lNg fc Romanil ([2004) find projected misalign- 
ments of ~ 1 0°, al thou gh most are consisten t with zero. 
iWang et all (|2006f ) and lNg fc Romanil (|2007f) find similar 
misalign ments for a larger sampl e using polarization data 
(also see iDeshpande et all 1 1 9991 ). but especially if we re- 
strict the sample to the younger pulsars (where Galactic 
acceleration should not have modifi ed the initial ve l ocity ) 
then the conclusions are simila r to lNg fc Romanil (2004) 
(also see Uohnston et aIll2005D . This suggests that, at 
least statistically, there still is a considerable degree of 
alignment between projected spin axes and proper mo- 
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tions. From this, IWang et all (|2007f ) and iNg fc Romanl 
argue that the asymmetries experienced by proto- 
neutron stars following core-collapse simultaneously can 
impart these stars with both kick and spin, and that 
these asymmetries consist of a stochastic ensemble of 
thrusts, each long enough to result in rotational averag- 
ing of the resultant linear momentum vector. The uncer- 
tainties discussed in this paper illustrate the difficulties 
in measuring precision alignments (or lack thereof) in 
any individual object. Further progress in these studies, 
for example, via detailed analyses of how the degree of 
alignment depends on parameters such as space velocity 
and surface magnetic field strength, is best achieved by 
adding to the total number of pulsars with information 
on the orientations of both spin and kick. 
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APPENDIX 
PROSPECTS FOR PARALLAX 

Given the importance of the Crab pulsar in our understanding of neutron stars, its precise distance remains a sur- 
prisingly open question. A trigonometric parallax has not yet been measured f or the pulsar. From th e dispersion of its 
radio pulses and a model of the Galactic electron density distribution ( NE2001.|Cordes fc L azio 2002), PSR B0531+21 
has an inferred distance of 1.7kpc (a distance range of 1.4 — 2.0kpc). iTrimbld (|1973f ) estimated a range of distances 
between 1.4 and 2.7 kpc based on a variety of lines of evidence, and the nominal distance to the Crab pulsar and 
its nebula is quoted as 2.0 ± 0.5 kpc. Precise measurements of the times of arrival of radio pulses have been used to 
measure radio pulsar parallaxes, but such measu rements require exceptio nal rotational stability, generally seen only 
in some recycled pulsars (e.g., PSR J 0437— 4715; Ivan Straten et aTll2001[ ). The young Crab pulsar has noisy timing 
residuals and shows rotational glitches (jWong. Backer, fc Lvnd l2001). ruling out such an approach to astrometry. Thus 
a parallax (and proper motion) for the Crab pulsar must rely on imaging at some wavelength range, and radio VLBI 
or optical observations with space telescopes are currently the most plausible approaches. 

From a purely numerical perspective, it appears that one should be able to use HST observations of the Crab 
pulsar, as described in this work, to measure its parallax. After all, f or bright stars the ePSF measurements and 
disto rtion solution are accurate to 0.01 pixel in an individual exposure (|Anderson fc Kindl2004 [2006; Ka plan et all 
12007ft . which is 0.25 mas for the ACS/High Resolution Camera (HRC) and 0.5 mas for the ACS/WFC. With its distance 
around 2 kpc, we expect a parallax near 0.5 mas, so with a sufficient number of exposures this should be measurable 
in principle. However, there are two limiting factors. First, it seems that for the brightest stars systematic effects 
prevent the combination of individual exposures fr om reducing the ast rometric uncertainty by the square-root of the 
number of exposures as one might expect (see, e.g.. lKaplan et al.ll2007l ). Second, unlike in radio interferometry where 
the parallax is measured relative to quasars and radio galaxies at essentially "infinite" distance, we must measure 
relative to other stars in our Galaxy that are at finite distances. For our measurement of the parallax of a neutron 
star at w 350 pc (jKaplan et al.l f2007h . we found that most of the reference stars were at 1-2 kpc and therefore had 
parallaxes at the 0.5-1.0 mas level. If we had ignored them, our parallax measurement would have been biased by the 
weighted mean of the parallaxes of the reference stars, or rs 0.5 mas, which is quite significant. Luckily, we had enough 
photometry of the field that we were able to determine photometric parallaxes for the reference stars. While not very 
accurate individually, they were sufficiently close to the true values (as we measured from our astrometry) to allow 
the correction of the ensemble of stars and the removal of the parallactic bias. 

From our color-magnitude diagram, we see that the background stars for the Crab pulsar are largely at 1-4 kpc, 
although there may be some at larger distances. Therefore the mean parallax of the background sources is likely 
> 0.25 mas, or approximately one half of the expected parallax of the Crab pulsar. This is a very significant correction 
and in order to measure the parallax of the pulsar with any significance we would need to know this bias to better 
than 10%. This might be possible with more detailed photometry (unfortunately, while the field is frequently observed 
by many facilities, most use narrow-band filters that are not suitable for spectral typing) or limited spectroscopy, but 
it would require a dedicated set of observations. 

Overall, then, the prospects for an optical astrometric parallax do not seem to be very good. First, we would need 
a number of additional astrometric HST observations, likely with the ACS/HRC (again, this may not be possible due 
to the recent failure of this instrument), where the Crab pulsar is not saturated, but this will of course mean that 
we detect fewer reference stars (due both to the limited field of the ACS/HRC and the shallower exposures). Since 
we are worried about accuracy at the < 0.01 pix level, we must also attempt to refine our estimate of the distortion 
solution and ePSF, which are difficult with relatively sparse fields like this one. We must also be confident in all of 
our systematics at this level. Second, we need a good number of multi-band photometric observations to measure 
reliable photometric parallaxes for at least the majority of the background stars; since we are observing out of the 
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Galaxy, we must be able to distinguish between solar metallicity st ars in the d i sk an d low-metallicity stars in the halo. 
The upcoming Wide Field Camera 3 (WFC3) on HST will have iStromgrenl (|1966l ) uvbyfi filters that greatly aid in 
stellar typing, but this may not be enough. It is therefore our opinion that an optical parallax is unlikely with current 
instruments, although it may be possible with future instruments such as the Space Interferometry Mission or GAIA. 

Very Long Baseline Interferometry (and specifically the Very Long Baseline Array) has been used to measure the 
proper motions and parallaxe s of a number of neutron stars, including s ome which are both weaker and more distant 
than the Crab pulsar (see, e.g. lBrisken et al.ll2002tlc"hatteriee et all 2005). Such observations would provide astrometry 
referenced to distant extragalactic quasars, eliminating uncertainty due to the reference frame (c.f. our discussion in 
§2.2, §2.5, and above), although we would still require DGR and LSR corrections. However, the Crab pulsar is 
embedded in an extremely radio-bright nebula, which dominates the system temperature of radio telescopes and 
thus limits the signal-to-noise ratio of radio interferometric observations. Coupled with the absence of a suitable 
extragalactic reference source nearby, the large increase in system temperature has limited attempts to measure a 
precise proper motion and parallax with the VLBA. Simply adding sensitivity (by increasing the collecting area, 
bandwidth, or integration time) does not address these limitations. We note that the ~6' size of the Crab nebula 
is comparable to the size of the primary beam of the 25-m VLBA antennas (~9' at 5 GHz). A telescope consisting 
of smaller dishes (such as the current Reference Design for the Square Kilometer Array 11 ) would have a wider field 
of view and suffer a proportionately smaller decrease in signal-to-noise ratio for the same total collecting area when 
observing a source as strong as the Crab nebula. The wider field of view and higher sensitivity would also provide 
suitable astrometric reference sources. Future radio telescopes with continent-sized baselines may thus enable a VLBI 
parallax for the Crab pulsar. 

It is unfortunate that the Crab pulsar, a subject of such intense and detailed investigation for so long, remains just 
beyond our current astrometric capabilities. However, the next generation of optical and radio telescopes should allow 
the measurement of a trigonometric parallax to this object, finally settling questions about its distance. 
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